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, The properties of spin polarized neutron matter are studied both at zero 

and finite temperature within the framework of the Brueckner-Hartree-Fock 
£^ ■ formalism, using the Argonne vl8 nucleon-nucleon interaction. The free en- 



ergy, energy and entropy per particle are calculated for several values of the 
spin polarization, densities and temperatures together with the magnetic sus- 
ceptibility of the system. The results show no indication of a ferromagnetic 
transition at any density and temperature. 
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Since the suggestion of Pacini [1] and Gold [2] pulsars are generally believed to be rapidly 
rotating neutron stars with strong surface magnetic fields in the range of 10 12 — 10 13 Gauss. 
Despite the great theoretical effort of the last forty years, there is still no general consensus 
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regarding the mechanism to generate such strong magnetic fields in a neutron star. The 
fields could be a fossil remnant from that of the progenitor star or, alternatively, they could 
be generated after the formation of the neutron star by some long-lived electric currents 
flowing in the highly conductive neutron star material. From the nuclear physics point of 
view, however, one of the most interesting and stimulating mechanisms which have been 
suggested is the possible existence of a phase transition to a ferromagnetic state at densities 
corresponding to the theoretically stable neutron stars and, therefore, of a ferromagnetic core 
in the liquid interior of such compact objects. Such a possibility has been considered since 
long ago by several authors within different theoretical approaches [3-23], but the results are 
still contradictory. Whereas some calculations, like for instance the ones based on Skyrme- 
like interactions predict the transition to occur at densities in the range (1 — 4)p (po — 0.16 
fm -3 ), others, like recent Monte Carlo [19] and Brueckner-Hartree-Fock calculations [20-22] 
using modern two- and three-body realistic interactions exclude such a transition, at least 
up to densities around five times p . This transition could have important consequences for 
the evolution of a protoneutron star, in particular for the spin correlations in the medium 
which do strongly affect the neutrino cross sections and the neutrino mean free path inside 
the star [24]. Therefore, drastically different scenarios for the evolution of protoneutron 
stars emerge depending on the existence of such a ferromagnetic transition. 

Most of the studies of the ferromagnetic transition in neutron and nuclear matter have 
been done at zero temperature. However, the description of protoneutron stars [25] motivates 
a study of spin polarized neutron matter at temperature T of the order of a few tens of MeV. 
Recently, the properties of polarized neutron matter both at finite and zero temperature, 
have been investigated [26] using a large sample of Skyrme-like interactions. The results 
of Ref. [26] indicate the occurrence of a ferromagnetic phase of neutron matter. However, 
contrary to what one would intuitively expect, the authors of Ref. [26] have found that 
the critical density at which ferromagnetism takes place decreases with temperature. This 
unexpected result was associated to an anomalous behaviour of the entropy of the system 
which becomes larger for the spin-polarized phase with respect the one for the non-polarized 
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phase, above a certain density. This was shown to be related to the dependence of the 
effective masses of neutrons with spin up and down on the amount of spin-polarization, 
and a new constraint on the parameters of the Skyrme force was derived if this anomalous 
behaviour is to be avoided [26] . 

In the present work, we study the bulk and single particle properties of spin-polarized neu- 
tron matter at finite temperature. To this aim we make use of a microscopic approach based 
on the Brueckner-Hartree-Fock (BHF) approximation of the Brueckner-Bethe-Goldstone 
(BBG) expansion. Here we make use of an extension of the BBG theory (i) to the case 
in which neutron matter is arbitrarily asymmetric in the spin degree of freedom [20] (i.e., 
Pl 7^ Pl, where pj (pjj is the density of neutron with spin up (down)), and (ii) to the 
case of finite temperature. In particular, we study the behaviour of the entropy of the sys- 
tem and the effective mass of neutrons as a function of the spin polarization parameter, 
A = (pj — Pi)/(p] + Pi)- We show that, contrary to what it is found in Ref. [26], the entropy 
of the polarized phase is lower than that of the non-polarized one, according to the idea that 
the polarized phase is more "ordered" than the non-polarized one. 

Our calculation starts with the construction of the neutron-neutron G-matrix, which 
describes in an effective way the interaction between two neutrons for each one of the spin 
combinations |1\ Tl, IT an d !!• This is formally obtained by solving the well known Bethe- 
Goldstone equation, written schematically as 

G (^)oi 0-2,0-30-4 — Kri02,0304 + J]] Koio 2 ,OiOj ' ! ^G ? (^)o i Oj,0304 5 (1) 

where the first (last) two sub-indices indicate the spin projection a =] (!) of the two 
neutrons in the initial (final) state, V is the bare nucleon-nucleon interaction, Q CiCj is the 
Pauli operator which allows only intermediate states compatible with the Pauli principle, 
and ui is the starting energy defined as the sum of the non-relativistic single-particle energies, 
£]{{)■> of the interacting neutrons. 

The single-particle energy of a neutron with momentum k and spin projection o —\ (!) 
is given by 
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s a (k) = y?- + Re[U a (k)} , (2) 

where the real part of the single-particle potential U a (k) represents the averaged field "felt" 
by the neutron due to its interaction with the other neutrons of the system. In the BHF 
approximation it is given by 

U a (k) = Y,n a ik')(kak'(x'\G{uj = e a (k) + e a ,(k')\kak'a') A , (3) 

a'k' 



where 



n a (k) = < 



1, if k < k a F 

(4) 

0, otherwise 



is the corresponding occupation number of a neutron with spin projection a and the matrix 
elements are properly anti-symmetrized. We note here that the so-called continuous prescrip- 
tion has been adopted for the single-particle potential when solving the Bethe-Goldstone 
equation. As shown by the authors of Refs. [28,29], the contribution to the energy per 
particle from three-body clusters is diminished in this prescription with respect to the one 
calculated with the gap choice for the single particle potential. We also note that the present 
calculation has been carried out using the Argonne vl8 nucleon-nucleon potential [30]. The 
momentum dependence of the single-particle spectrum can be characterized by the effective 
mass m* a (k) defined as: 

ml(k) _ k f deJJ^ y 1 
m m\ dk J 

where m is the bare neutron mass. 

The total energy per particle is easily obtained once a self-consistent solution of Eqs. 
(l)-(3) is achieved 

The many-body problem at finite temperature has been considered by several authors 
within different approaches, such as the finite temperature Green's function method [31], 



the thermo field method [32], or the Bloch-De Domicis (BD) diagrammatic expansion [33]. 
The latter, developed soon after the Brueckner theory, represents the "natural" extension 
to finite temperature of the BBG expansion, to which it leads in the zero temperature limit. 
Baldo and Ferreira [34] showed that the dominant terms in the BD expansion were those 
that correspond to the zero temperature of the BBG diagrams, where the temperature is 
introduced only through the Fermi-Dirac distribution 

UK T) = l + exp([e a (k,T)-^(T)]/T) ' (?) 

HaiT) being the chemical potential of a neutron with spin projection a. Therefore, at the 
BHF level, finite temperature effects can be introduced in a very good approximation just 
replacing in the Bethe-Goldstone equation: (i) the zero temperature Pauli operator Q aiCTj = 
(1 - n (T J(l - n aj ) by the corresponding finite temperature one Q aiaj (T) = {1- f a .)(l- f a .), 
and (ii) the single-particle energies e a (k) by the temperature dependent ones e a (k,T) ob- 
tained from Eqs. (2) and (3) when n a (k) is replaced by f a (k,T). These approximations, 
which are supposed to be valid in the range of densities and temperatures considered here, 
correspond to the "naive" finite temperature Brueckner-Bethe-Goldstone (NTBBG) expan- 
sion discussed in Ref. [34]. 

In this case, however, the self-consistent process implies that, together with the Bethe- 
Goldstone equation and the single-particle potential, the chemical potentials of neutrons 
with spin up and down must be extracted at each step of the iterative process from the 
normalization condition 

P, = E/^ T )' (8) 

k 

This is an implicit equation which can be solved numerically. Note that the G-matrix ob- 
tained from the Bethe-Goldstone equation (1) and also the single-particle potentials depend 
implicitly on the chemical potentials. 

Once a self-consistent solution is achieved the total free energy per particle is determined 

by 



where E/A is evaluated from Eq. (6) replacing n a (k) by f a (k,T) and the total entropy per 
particle, S/A, is calculated through the expression 

| = -\T,[Uk,T)\n(f a (k,T)) + (1 - f a (k,T)Ml - MKT))] . (10) 

From the free energy per particle, we can get the remaining macroscopic properties of the 
system. In our case, we are particularly interested in the magnetic susceptibility Xi which 
characterizes the response of a system to a magnetic field and gives a measure of the energy 
required to produce a net spin alignment in the direction of the field. It is given by 

Y = £i (id 

X ( d 2 {F/A) \ ^ LL > 

where fi is the magnetic moment of the neutron. 

The single-particle potentials of neutrons with spin up and down have been simultane- 
ously and self-consistently calculated together with their effective interactions. The results 
at p = 0.16 fm -3 and spin polarization A = 0.5 are reported for T=0 (left panel) and T=40 
MeV (right panel) on the top panels of Fig. 1. The neutron single-particle potential splits 
up in two different components when a partial spin polarization is assumed. In the case 
of Fig. 1, the single-particle potential Re[U^(k)} for neutrons with spin up (the most abun- 
dant component) is less attractive than the one for neutrons with spin down, Re[U±(k)]. 
As demonstrated by the authors of Ref. [21] (see in particular their Eqs. (23) and (24)), 
this splitting (i) is the result of a phase space effect, i.e. to the change in the number of 
pairs which the neutron under consideration \k,a) can form with the remaining neutrons 
\k < hp, a' =t, I) of the system as neutron matter is polarized, and (ii) is due to the spin 
dependence of the neutron-neutron G-matrix in the spin polarized medium (see Eq. (1)). 
Indeed, as polarization increases, the single particle potential of a spin up neutron is built 
from a larger number of up-up pairs that form a spin triplet state (S = 1) and, due to 
the Pauli principle, can only interact through odd angular momentum partial waves. Con- 
versely, the potential of the less abundant species is built from a relatively larger number 
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of up-down pairs which can interact both in the S = and S = 1 two body states. Thus, 
the potential of the less abundant species receives also contributions from some important 
attractive channels as e.g. the 1 S'o. 

The increase of the temperature changes moderately the single-particle potentials. The 
real part becomes slightly less attractive, whereas the imaginary part increases in size as a 
consequence of the increase of phase space in the low momentum region. 

The momentum dependence of the corresponding effective masses of the two components 
is also shown in the bottom panels of the figure for the same values of density, spin polariza- 
tion and temperatures. The general effect of temperature is to smooth out the enhancement 
of the effective mass near the Fermi surface, as observed in the work of Ref. [35] in symmetric 
nuclear matter. 

In Fig. 2 we show the effective mass m*(m|) for neutron with spin up (down) as a 
function of the spin polarization A, for fixed density (p = 0.16 fm~ 3 ) and temperature (T=0 
and T=40 MeV). The effective mass is calculated using Eq. (5) taken for each component 
at the corresponding Fermi momentum. Obviously, for A = the effective mass of the 
two components coincides. Once some amount of polarization is considered, the values of 
the effective masses split in two, the effective mass of the most abundant component being 
larger than the one of the less abundant. As can be seen the effective masses show an almost 
linear and symmetric variation with respect to their common value at spin polarization 
A = 0, both at T=0 and T=40 MeV. Deviations from this behaviour are only found at the 
higher polarization values. This behaviour of m* is a direct consequence of the scissors- 
like dependence of the single particle potential Re[U a ] as a function of the spin polarization 
parameter A (see Fig. 2 of Ref. [21]). A similar qualitative behaviour for the nucleon effective 
mass, as a function of the isospin asymmetry parameter, (3 = (p n — p p ) / p, has been found 
in isospin asymmetric nuclear matter [36-38] (see in particular Eq. (94) in Ref. [37]). 

The differences of the free energy (F/A), energy (E/A) and entropy (S/A) per particle 
between the totally polarized and the non-polarized phases are reported in the left, central 
and right panels of Fig. 3 as a function of the density for several temperatures. The differ- 
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ences in the three quantities increase with density and increase (decrease) with temperature 
in the case of the free energy (energy and entropy). Contrary to the results of Ref. [26] 
with the Skyrme interaction, these differences are always positive for the F/A and E/A. 
This is an indication that the non-polarized phase is energetically preferred in the range of 
densities explored. Therefore, we can conclude that a phase transition to a ferromagnetic 
state is not to be expected from our microscopic calculation. If such a transition would 
exist, the difference in the free energy would become zero at some density, indicating that 
the ground state of the system would be ferromagnetic from that density on. In addition, 
the difference in the entropy is always negative indicating, as one intuitively expects, that 
the totally polarized phase is more "ordered" than the non-polarized one. 

In Fig. 4 we show the behaviour of the free energy F/A per particle as a function of 
the spin polarization for several densities (left panel) and temperatures (right panel). Cir- 
cles, squares, diamonds and triangles correspond to our BHF results, whereas the solid lines 
correspond to the parabolic approximation discussed below. As we expected from our pre- 
vious calculations at zero temperature [20] and [21], F/A is symmetric in A and it shows 
a minimum at A = for all the densities and temperatures considered. This is again an 
indication that the ground state of neutron matter is paramagnetic, in opposition to what it 
is found in Ref. [26] for Skyrme-like interactions where, as a consequence of the anomalous 
behaviour of the entropy, the minimum of F/A is situated at < A < 1 and moves to 
higher polarizations when the temperature increases. It is also interesting to note that the 
dependence of F/A on the spin polarization is "up to a very good approximation" parabolic. 
One can try to characterize that dependence in the following simple analytic form: 

^(p,A,T) = ^(p,0,T) + a(p,T)A 2 (12) 

where, assuming the quadratic dependence to be valid up to |A| = 1 as our results indicate, 
the value of a(p, T) can be easily obtained for each density and temperature as the difference 
between the total free energies per particle of totally polarized and non-polarized neutron 
matter 
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a(p,T) = ^(p,±l,T)-j(p,0,T) . (13) 

The magnetic susceptibility can be evaluated then in a very simple way if the parabolic 
dependence of Eq. (12) is assumed, giving 

^• T »-^T)- < 14 > 

The ratio Xf/Xi where xf is the magnetic susceptibility of the free Fermi gas, is shown in 
Fig. 5 as a function of density for several temperatures. Starting from 1, the ratio increases 
as the density increases at any temperature and no signal of a change of such a trend is 
expected at higher densities, contrary to the results of Ref. [26] in the case of the Skyrme- 
like interactions. This is again an indication that a ferromagnetic transition, whose onset 
would be signaled by the density at which this ratio becomes zero, is not seen and not 
expected at larger densities either. 

Finally, the behaviour of the entropy per particle S/A as a function of the spin polar- 
ization at a fixed density p = 0.32 fm~ 3 for several temperatures is shown in Fig. 6. The 
entropy, as the free energy, is also symmetric and almost parabolic in A. Its maximum 
is placed at A = for all the densities and temperatures considered, as one naively ex- 
pects, contrary to the findings of Ref. [26]. In this reference, it was shown that for a pure 
parabolic single particle spectrum, as it is the case for the Skyrme interaction, imposing 
the entropy of the polarized phase to be smaller than the unpolarized one for a given den- 
sity and temperature, is equivalent to requiring the ratio of the neutron effective masses in 
the fully polarized and unpolarized phases to be smaller than 2 2 / 3 . In the BHF approach, 
the momentum and temperature dependence of the effective mass prevents from deriving a 
similar rigorous condition. However, thinking in terms of a value of the effective mass that 
would characterize the single particle spectrum in average, or considering just the effective 
mass at the Fermi surface, which is the most relevant for the calculation of the entropy 
at small temperatures, we can then explore if the BHF calculations respect the condition 
derived in [26]. In fact, in the case of p = 0.16 fm~ 3 and T = 40 MeV we find (see Fig. 2) 
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m*(A = l)/mfQ)(A = 0) = 1.09, which is smaller than the limit established in Ref. [26]. 
This is true for all the densities and temperatures explored in this work and therefore the 
entropy of the polarized phase is always smaller than that for the unpolarized one. 

In summary, we have studied the properties of spin polarized neutron matter both at 
zero and finite temperature within the framework of the Brueckner-Hartree-Fock formalism. 
We have determined the single-particle potentials and the effective mass of neutrons with 
spin up and down for arbitrary values of the density, temperature and spin polarization. 
We have found that the spin up and spin down effective masses show an almost linear and 
symmetric variation with respect to their values at spin polarization A = 0. 

We have determined the differences of the free energy (F/A), energy (E/A) and entropy 
(S/A) per particle between the totally polarized and non-polarized phases. We have found 
that, in contrast to the results of a similar study with the Skyrme interaction [26], these 
differences are always positive for the F/A and E/A which is an indication that the non- 
polarized phase is energetically favorable, from which we can conclude that a phase transition 
to a ferromagnetic state is not to be expected. In addition, contrary to the results with the 
Skyrme interaction, we have found that the difference in the entropy is always negative 
according to the idea that the totally polarized phase is more "ordered" than the non- 
polarized one. 

Finally, we have seen that both the free energy and the entropy per particle are not only 
symmetric on the spin polarization but also parabolic in a very good approximation up to 
|A| = 1. This finding supports the calculation of the magnetic susceptibility by using only 
the free energies of the fully polarized and non-polarized phases. 

This work is partially supported by DGICYT (Spain) project BFM2002-01868 and by 
the Generalitat de Catalunya project 2001SGR00064. This research is part of the EU In- 
tegrated Infrastructure Initiative Hadron Physics project under contract number RII3-CT- 
2004-506078. One of the authors (A. Rios) acknowledges the support from DURSI and the 
European Social Funds. 



10 



REFERENCES 

[1] F. Pacini, Nature (London) 216, 567 (1967). 
[2] T. Gold, Nature (London) 218, 731 (1968). 

[3] D. H. Brownell and J. Callaway, Nuovo Cimento 60 B, 169 (1969). 
[4] M. J. Rice, Phys. Lett. 29 A, 637 (1969). 

[5] J. W. Clark and N. C. Chao, Lettere Nuovo Cimento 2, 185 (1969). 

[6] J. W. Clark, Phys. Rev. Lett. 23, 1463 (1969). 

[7] S. D. Silverstein, Phys. Rev. Lett. 23, 139 (1969). 

[8] E. Cstgaard, Nucl. Phys. A 154, 202 (1970). 

[9] J. M. Pearson and G. Saunier, Phys. Rev. Lett. 24, 325 (1970). 
[10] V. R. Pandharipande, V. K. Garde and J. K. Srivastava, Phys. Lett. 38 B, 485 (1972). 
[11] S. O. Backman and C. G. Kallman, Phys. Lett. 43 B, 263 (1973). 
[12] P. Haensel, Phys. Rev. C 11, 1822 (1975). 

[13] A. D. Jackson, E. Krotscheck, D. E. Meltzer and R. A. Smith, Nucl. Phys. A 386, 125 
(1982). 

[14] M. Kutschera and W. Wojcik, Phys. Lett. B 223, 11 (1989). 

[15] S. Marcos, R. Niembro, M. L. Quelle and J. Navarro, Phys. Lett. B 271, 277 (1991). 
[16] P. Bernardos, S. Marcos, R. Niembro, M. L. Quelle, Phys. Lett, B 356, 175 (1995). 
[17] A. Vidaurre, J. Navarro and J. Bernabeu, Astron. Astrophys. 135, 361 (1984). 
[18] M. Kutschera and W. Wojcik, Phys. Lett. B 325, 271 (1994). 
[19] S. Fantoni, A. Sarsa and K. E. Schmidt, Phys. Rev. Lett. 87, 181101 (2001). 

11 



[20] I. Vidaha, A. Polls and A. Ramos, Phys. Rev. C 65, 035804 (2002). 
[21] I. Vidaha and I. Bombaci, Phys. Rev. C 66, 045801 (2002). 

[22] W. Zuo, U. Lombardo and C.W. Shen, in Quark-Gluon Plasma and Heavy Ion Collisions, 
Ed. W.M. Alberico, M. Nardi and M.P. Lombardo, World Scientific, p. 192 (2002). 

[23] A. A. Isayev and J. Yang, Phys. Rev. C 69, 025801 (2004). 

[24] J. Navarro, E. S. Hernandez and D. Vautherin, Phys. Rev. C 60, 045801 (1999). 

[25] M. Prakash, I. Bombaci, Manju Prakash, P.J. Ellis, J.M. Lattimer and R. Knorren, 
Phys. Rep. 280, 1 (1997). 

[26] A. Rios, A. Polls and I. Vidaha, Phys. Rev. C 71, 055802 (2005). 

[27] J. Margueron, J. Navarro and N. V. Giai, Phys. Rev. C 66, 014303 (2002). 

[28] H. Q. Song, M. Baldo, G. Giansiracusa and U. Lombardo, Phys. Rev. Lett. 81, 1584 
(1998). 

[29] M. Baldo, G. Giansiracusa, U. Lombardo and H. Q. Song, Phys. Lett. B 473, 1 (2000). 

[30] R. B. Wiringa, V. G. J. Stoks, R. Schiavilla, Phys Rev. C 51, 38 (1995). 

[31] A. L. Fetter and J. D. Walecka, Quantum Theory of Many Particle Physics (McGraw- 
Hill, New York, 1971). 

[32] For a recent review see P. A. Hening, Phys. Rep. 253, 235 (1995). 

[33] C. Bloch, Nucl. Phys. 7, 451 (1958); C. Bloch and C. De Domicis, ibid 7, 459 (1958); 
10, 181 (1959); 10, 509 (1959). 

[34] M. Baldo and L. Ferreira, Phys. Rev. C 59, 682 (1999). 

[35] A. Lejeune, P. Grange, M. Martzolff and J. Cugnon, Nucl. Phys. A 453, 189 (1986). 
[36] I. Bombaci and U. Lombardo, Phys. Rev. C 44, 1892 (1991). 

12 



[37] I. Bombaci, in Nuclear methods and the nuclear equation of state, Ed. M Baldo, World 
Scientific, p. 381 (1999). 

[38] B.-A. Li, Phys. Rev. C 69, 064602 (2004). 



13 



120 



> 

J3 40 





D 

-40 
-80 
1 

0.95 
0.9 



1 1 1 


| 1 | 1 | 1 | 1 | 1 


- T=0MeV,A=0.5 Re /- 




/ / 


■ I I 










— - 






1 . 1 . 1 . 1 . 1 . " 



o 0.85 
_> 

0.8 

pq 

0.75 



I 1 I 1 I 1 I 1 I 




.4 

l ■ l ■ l 



- neutron up 
■ neutron down 
I ■ I ■ I ■ 



k [fin ] 



I 1 I 



T 



4 5 
I 1 I 



6 

T 



7 



T=40 MeV, A=0.5 

-3 




I ■ I ■ I ■ I ■ I ■ I 




120 
80 
40 


-40 
-80 
1 

0.95 

0.9 

0.85 

0.8 

0.75 



k[fm ] 



FIG. 1. Single-particle potential (top panels) and effective mass (bottom panels) of neutrons 
with spin up (solid lines) and spin down (dashed lines) as functions of the linear momentum at 
fixed density (p = 0.16 fm~ 3 ) and spin polarization (A = 0.5) for T=0 (left panels) and T=40 
MeV (right panels). The arrows denote the value of the corresponding Fermi momenta. 
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FIG. 2. Neutron effective mass at the corresponding Fermi surface of the spin up and down 
components as a function of the spin polarization at p = 0.16 fm -3 for T = and T = 40 MeV. 
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FIG. 3. Differences of the free energy (left panel), energy (central panel) and entropy per 
particle (right panel) between fully polarized and non-polarized neutron matter as a function of 
density for several temperatures. 
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FIG. 4. Left panel: free energy per particle at zero temperature as a function of the spin 
polarization for several densities. Right panel: free energy per particle at a fixed density p = 0.36 
fm~ 3 as a function of the spin polarization for several temperatures. Circles, squares, diamonds and 
triangles show our BHF results, whereas the solid lines correspond to the parabolic approximation 
defined in Eq. (12). 
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FIG. 5. Ratio between the magnetic susceptibility of the free Fermi gas and the corresponding 
magnetic susceptibility of interacting neutron matter as a function of density for several tempera- 
tures. 
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FIG. 6. Entropy per particle as a function of the spin polarization at p = 0.32 fm 3 for several 
temperatures. 
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